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ABSTRACT 

We report a scaling of the mean matter density with the width of the saturated 
Lya absorptions. This property is established using the "pseudo-hydro" technique 
ijCroft et alJ ll9981. It provides a constraint for the inversion of the Lya forest, which 
encounters difficulty in the saturated region. With a Gaussian density profile and the 
scaling relation, a simple inversion of the simulated Lya forests shows that the one- 
dimensional mass power spectrum is well recovered on scales above 2 ft, Mpc, or 
roughly k < 0.03 km _1 s, at z = 3. The recovery underestimates the power on small 
scales, but improvement is possible with a more sophisticated algorithm. 
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1 INTRODUCTION 

The Lya forest has offered a unique way to study the large- 
scale structure of the universe at high redshifts. It is based on 
the theory that the Lya forest is a result of the absorptions 
of quasar continuum by the diffusely distributed and pho- 
toionized intergalactic medium (IGM), and the IGM traces 
the density of the und erlying mass field on scales larger than 
the Jeans length llBlJTiffli ICen et al.lll994l : iBi et al.lll99l 
iPetitiean et alJll99St Iviel et al]l20o3h In other words, the 
baryon density in Lya absorption systems is roughly propor- 
tional to the matter dens i ty, when p < 10 |Bi fc Davidsenl 
ll997l:lGnedin fc Huill998tlZhang et allll998ft .~ 

To extract the baryon density pb or the mass density 
p, one may invert the transmitted flux F of the Lya forest 
using 



F 



-A* 



-Apt 3 



(1) 



where: r is the Lya optical depth; A depends on the tem- 
perat ure of IGM, the Hu photoionization rate, and red- 
shift jRauch et alJll997t ICroft et al.lll99Sft: and 3. rela ted 
to equation of state of the IGM iHui fc Gnedinlll997lh is 
in the range 1.6-1.8 JCroft et al.lll998ft . The densities are in 
units of their corresponding cosmic mean density, and so the 
overdensity S = p — 1 . 

When the density is high enough, the spectrum is sat- 
urated, i.e. F ~ 0. With noises and uncertainties in the 
spectrum, the direct inversion using equation Q is very un- 
reliable in the saturated region. Despite the difficulty, meth- 
ods of direct inversion are systematically developed, for ex- 
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ample, with Lucy's method bv iNusser fc HaehnebJ lll999ft . 
and wi th Bayesian method for a three-dimensional inver- 
sion bv lPichon et alJ (1200 if) . One may also use higher order 
lines to recover the optical depth and the unde rlying density 
JCowie fc Songaildll99dlAguirre et al.ll 2002ft . even though 
the contamination by lower order lines needs to be carefully 
removed. Once the cosmic density field is obtained, many 
statistics, such as the mass power spectrum, can be mea- 
sured. 

A nother approach is taken bv lCroft et alJ <ll998t IT9991 . 

2002) to map the mass power spectrum directly from the 
flux power spectrum of the Lya forest without an inversion. 
Thus the saturation problem is avoided. However, a close 
examination of the Fourier transform of equation Q shows 
that powers on different scales are mixed by the non-linear 
density-flux relation (see Section • The mixing depends 
on the underlying density field, and it is hard to predict 
analytically. 

If the inversion is necessary, a proper treatment in the 
saturated region has to be developed. In many physical sys- 
tems, sizes are often correlated with other quantities such 
as masses and densities. For example, more massive stars or 
dark matter haloes have larger sizes, but lower mean den- 
sities jBinnev fc Merrifieldll99l iNavarro et all 1996ft . One 
may expect a similar trend for the saturated Lya absorp- 
tion. On the contrary, the mean density is found to increase 
with the width of saturation. This is due to the fact that 
the IGM is very diffuse and far from virialization, while the 
other objects mentioned are the opposite. 

This paper is organised as follows. Section presents 
the scaling relation. A simple inversion with the scaling is 
tested in Section|3for both the one-dimensional mass power 
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spectrum and the one-point distribution function. Section 
0| examines the mapping between the mass power spectrum 
and the flux power spectrum. The discussion and conclusions 
are given in Section^] Unless addressed otherwise, the power 
spectra below are implicitly one-dimensional. 



2 THE SCALING 

Neglecting the probability that two physically separate ab- 
sorption systems fall in the same redshift coordinate, one 
can associate a saturated absorption in the Lya forest with 
a single high density region. If we assume further that there 
is no substructure present, and the IGM evolves more or 
less the same way everywhere, then the size of the saturated 
region has to be tightly correlated with its mean density 
<Schavdl200l . In reality, the neglected elements above and 
uncertainties elsewhere will introduce a spread to the corre- 
lation. 



2.1 Simulated Lya Forest 

It is demonstrated by various hydrodynamic a l simulations 



Table 1. Parameters of the N-body simulations. 



it is demonstrated by various hydrodynamical simulations 
iBi fc DavidserJ Il997t iGnedin fc Hu] Il998t IZhang et alJ 



119981: iBrvan et alJll99a iDave et alJll999t l200lh that the 
baryon density, the H I column density, and the IGM temper- 
ature correlate well with the underlying dark matter density 
for p < 10. This leads to the useful "pseudo-hydro" tech- 
nique JCroft et alJll998t) for simulating the Lya forest with 
collisionless N-body simulations instead of hydrodynamical 
simulations. Basically it draws the line-of-sight (LOS) dark 
matter density in redshift space from N-body simulations, 
and then converts the density to flux using equation (0. 

One may have concerns about the Jeans length, below 
which the pseudo-hydro technique may not be applicable. 
The comoving Jeans length Lj is 



87T7fcBT 



s(l+f) 

VGp 



3pm p H^Qp(l + z) 



440 /i _1 kpc T 4 1/2 [f7p(l + z)] 



-1/2 



(2) 



where c s is the speed of sound, G is the gravitational con- 
stant, 7 = 5/3 is the ratio of specific heats, ks is the 
Boltzman constant, T is the temperature, n — 0.59 is the 
mean molecular weight, m p is the proton mass, Ho — 100 
km s _1 Mpc -1 , fi is the matter density in units of the crit- 
ical density at present, and T4 ~ 1.5 is the temperature in 
units of 10 K. It is understood that the first p in equa- 
tion is the actual density, while the rest are densities 
in units of the cosmic mean density. For the low-density 
cold dark matter (LCDM) model at z — 3 and p = 1, 
L.j = 490 h~ 1 kpc. The Jeans length suggests that the small- 
est scale we can analyse reliably without hydrodynamical 
simulations is 0.5 /i -1 Mpc. However, saturated regions usu- 
ally have 10 > p > 5 (see Fig.0, and the tem perature scales 
with the density as T oc p a with a = 0.3-0.6 llHui fc Gnedinl 
119971) . so the Jeans length i s reduced by a factor o f 1.4-2.2 in 
those regions. In addition, IGnedin fc Huil l[l998) point out 
that the actual linear filtering scale is 1.5-2.5 times smaller 
than the Jeans length at z = 3, such that baryons trace the 
dark matter on scales even half of Lj. Hence, the pseudo- 
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h is the Hubble constant in units of 100 km s _1 Mpc -1 . 

r is the shape parameter in the power spectrum. 

<T8 is the r.m.s. density fluctuation within a radius of 8 /i —1 Mpc. 

hydro technique is sufficient for the analysis above the scale 
of 200 /i _1 kpc. 

A standard particl e-particle-particle- mesh (P 3 M, 
iHocknev fc Eastwoodlll98lh code developed bv l.Iing fc Fanel 
( 1994) is used to evolve 128 3 dark matter particles in a 
cubic box of 12.8 ft _1 Mpc (comoving) each side. The 
initi al power s p ectru m is given by the fitting formula from 
iBardeen et alJ l)l98rtl . The model parameters are listed in 
Table All the models start from z — 15, and stop at 
2 = 3 in 950 steps. 

The parameter A in equation Q is chosen to fit the 
mean flux (F(z)) = exp[-r(z)1, wh ere t(z) = 0.0032(1 + 
z) 3 - 37 for z = 3-4 llKim et alJl2002l) . This mean flux for- 
mula is consistent with other observation s l)Lu et alJll99rl : 
iRauch et al"]ll997l : iMcDonald et all 120001) . We extrapolate 
the mean flux up to 2 = 4.5, which is not critical to the 
analysis i n Secti on l2.2l bu t is n evertheless supported by sim- 
ulations (Riedi eer et al]|l998T) . The constant /3 in equation 
Q is set to 1.6. Four samples of the Lya forest are shown in 
Fig. IH They are drawn from the LCDM simulation at 2 = 3. 
The corresponding LOS densities are plotted along with the 
fluxes. It is evident that most of the Lya lines arise where 

10. 



2.2 The Scaling 

We define the width w of the saturation as the distance 
between the two points that bracket the absorption at a 
given threshold flux level r\. In other words, it is the width 
of a region in which F < rj. Since the flux F in the saturated 
region is dominated by noise, one should set the threshold 
above the noise level. Within a reasonable range of spectral 
quality, we choose the threshold r\ — 0.03 and 0.10 for the 
tests below. 

Once the width is determined, the mean density p can 
be readily calculated from the corresponding LOS density 
field. Since the density and the flux are assigned only on 
grid points, interpolation is needed to find w and p. Fig. [5] 
shows the scaling relation for the LCDM model at 2 = 3. It 
clearly demonstrates that p increases with w. Furthermore, 
the correlation is reasonably tight for narrow saturations. 
Notice that the trend lines are obtained by fitting only the 
data with w ^ 2 h~ Mpc, because wider saturations are 
very rare as compared to the rest. 

The p-w relation is well fitted by 



p = po + a (w/h """Mpc) 6 



(3) 



The term po sets a baseline for the scaling, because the mat- 
ter density is non-vanishing even if there is no saturation. 
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Figure 1. Samples of the simulated Lya forest from the LCDM 
model at z = 3. The solid lines are the flux F, the dashed lines are 
the matter density p, the dotted lines are the recovered density, 
and the horizontal dash-dotted lines indicate the threshold flux 
level r\ = 0.03 (see Section[3J in the upper 2 panels, and r\ = 0.10 
in the lower 2 panels. The scale of the densities, in units of cosmic 
mean density, is indicated on the right axis. The flux and the 
density are both presented in redshift space. 



Table 2. Parameters in the fitting p = po + a (w/h ^-Mpc) 6 for 
the LCDM model. 



z 


'/ 


PO 


a 


b 


3 


0.03 


1.93 


3.85 


0.545 




0.10 


1.39 


3.39 


0.564 


4.5 


0.03 


0.11 


2.78 


0.455 




0.10 


0.24 


2.04 


0.556 



The parameter a is a scaling factor, which reflects the over- 
all amplitude of the density fluctuation, and the exponent 
b is more or less determined by the nature of hierarchical 
clustering. The width w is in units of /i _1 Mpc. Table [5] lists 
the values for LCDM model at z = 3 and 4.5. The reason 
why po decreases with redshift is that the universe is more 
uniform early on, so that even low density regions have to 
absorb a substantial amount of Lya flux to produce the low 
mean flux. This is possible because the neutral fraction of 
the IGM at z — 4.5 is higher than that at z = 3. As the 
universe evolves, the density fluctuation grows stronger and 
stronger, and the scaling factor a becomes larger and larger. 
The exponent b has changed little over z = 3-4.5. The pa- 
rameters also show a dependence on the threshold flux r/, 
because r\ sets the threshold density p v , above which the 
p-w relation is explored. 
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Figure 2. The correlation between the mean density p and the 
width w of the saturated region in the Lya forest. The trend lines 
fit the data well beyond 2 h~ 1 Mpc, even though they are obtained 
by fitting only the data with w $J 2 h~ 1 Mpc. The data above 
0.3 /i- 1 Mpc are binned in logarithmic intervals for readability. 
Similar treatment applies to Figs |1] |S| and 151 as well. 
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Figure 3. The best-fittings of the p-w relation for the 4 models at 
z = 3-4.5. The threshold r\ = 0.03. The parameters are obtained 
by fitting the saturations with w ^ 2 h~ 1 Mpc. 

Fig. |S| shows the evolution of the p-w relation for the 
4 models. The difference between models is mostly due to 
the amplitude of fluctuations - in other words, the power 
spectrum. 

2.3 The Physics 

The p-w relation is analogous - bu t not completely equiv - 
alent - to the curve of growth (e.g. iPress fc Rvbickilll993Tl . 
which studies the correlation between the H I column density 
iVm and the equivalent width W of Lya absorptions. The 
similarity is as follows. The width w is approximately the 
same as W, and the mean matter density p is proportional 
to to _1 7Vhi if p oc pb holds true. Since iVm oc W at small 
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values of W, it is not surprising to see p grow slowly, i.e. 
dlnp/dlnw ~ 0, for small w. 

On the other hand, the p-w relation addresses saturated 
absorptions where, according to the curve-of-growth analy- 
sis, W is almost a constant independent of the Hi column 
density. Therefore, p would have risen steeply against w, i.e. 
b 3> 1. This apparent inconsistency arises from the cosmo- 
logical context, because the width w is determined not only 
by Nm but also by the physical extent of the (relatively) 
dense region through Hubble expansion and peculiar veloc- 
ities. 

Although it is not obvious why b is less than unity, we 
can make an order-of-magnitude estimate using the density 
profile p(r) oc r _1 from a spherical self-similar infall. We 
modify the profile to avoid the singularity at r = by adding 
a smoothing length e, so that p(r) oc (r 2 + e 2 ) - • The mean 
density within the radius R, at which p = p v , is 



0.01 



0.1 (km" 1 s) 



p(r)r dr = 



3p v VR 2 + 



27? 3 



(4) 



R\/R 2 + e 2 + e 2 log 



From equation (QJ, one gets an estimate b ~ dlnp/dlni? = 
0.39, 0.42, and 0.43 for e = 0.2, 0.3, and 0.4 ft _1 Mpc re- 
spectively at R — 2 /i _1 Mpc. The same quantity for the 
p-w relation is din p/d\nw — 0.41, where we have used the 
parameters of the LCDM model at z = 3 with r\ = 0.03 
(see Table |5J. It should be stressed that the p(r) oc r~ 
profile is not quite justified for the IGM at z = 3, and with- 
out any modification, it gives the same mean density within 
the boundary p — p v regardless the size of the system, i.e. 
b = 0. Therefore, equation @ is not expected to give a fair 
approximation of the p-w relation. 



3 INVERSION WITH A GAUSSIAN DENSITY 
PROFILE 

If the flux level is above the threshold, equation Q can be 
used to find the density. When it is below the threshold, one 
must provide a density profile that matches p(w) and p n to 
fill in the missing information in the saturation. 

While the density profile is worth studying in its own 
right, we simply choose the Gaussian profile 

B 



P(s) 



2tt 



cxp 



1 



(5) 



where s is the coordinate in redshift space, so is the centre 
of the saturation, and B and a are solved simultaneously 
from 



p(w) 



u/2 



u/2 



p(s)ds = — erf(^-= — Y 
w V2V2 a) 



(6) 



Pv = p(s ± w/2) = 



B 



exp 



1 /W\2 

"8V' 



where erf(a;) is the error function. The Gaussian profile 
has the advantage that it does not introduce any artificial 
power on small scales. However, it is arguable that the right 
amount of small-scale power should be added through the 
profile, and so a more realistic profile may be needed. 

The recovered LOS densities with 77 = 0.03 and r\ — 0.10 
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Figure 4. The power spectra of original densities and recovered 
densities in redshift space. Circles include all LOS densities from 
the N-body simulation, and the dashed line contains only the 
ones that have a maximum width of saturation wm < 1 1 Mpc. 
The recovered densities are inverted from fluxes with thresholds 
r\ = 0.03 (solid line) and r\ = 0.1 (dotted line). The error bars 
are la dispersions among 610 groups, each of which consists of 
20 LOS densities. The error bars of the recovered densities, which 
are not plotted, are comparable to that of the original densities. 



are shown along with the four original LOS densities and 
fluxes in Fig. Q Since a universal density profile is employed 
in the inversion, the recovered densities do not necessarily 
match the original densities. 

To assess the statistical quality of the inversion, we plot 
in Fig. 0]thc original and the recovered mass power spectra 
of the LCDM model with different flux thresholds. Other 
models are shown in Fig. The recovered power spectrum 
agrees well with the original power spectrum on large scales 
(k < 3 h Mpc -1 ), but it is underestimated on smaller scales, 
where the Gaussian profile essentially has no power. The 
signal to noise ratio, or the threshold flux, has little influence 
on large scales, but a low noise level does slightly improve 
the recovery on small scales. 

For comparison, we show in Figs 2] and |^| the original 
power spectrum of the LOS densities that have a maximum 
width of saturation wm < 1 /i~ 1 Mpc. This is equivalent to 
removing - or giving less weight to - the saturated regions 
when one measures the power spectrum. It is seen that the 
saturated regions are very important to scales > 2 ft _1 Mpc, 
while the unsaturated regions give a good estimate of the 
small-scale power. Thus one may improve the recovery of 
power spectrum as follows. First, invert the Lycv forest with 
equation Q and the p-w relation. Second, do the inver- 
sion after removing the saturated regions. Finally, the best- 
estimate of the power spectrum is just the common envelope 
of the power spectra of the two recovered densities. 

Fig. [^indicates that the p-w relation varies from model 
to model. Thus it is necessary to check if the recovery is sen- 
sitive to p(w) - in other words, if it is model dependent. In 
Fig. [(Jwe plot the power spectra of LCDM densities recov- 
ered with po(w) and Pt(w) from open cold dark matter 
(OCDM) and tilted cold dark matter (TCDM) models re- 
spectively. It seems that by boosting p(w) a small amount 
[Pt(w) lies a little higher than pl(w)], one gets even better 
estimate of the power spectrum on small scales. However, 
lowering p(w) could underestimate the power by a factor of 
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Figure 5. The ame as Fig.H but for OCDM, SCDM, and TCDM 
models. 
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Figure 6. The same as Fig. |1] except that the fluxes are drawn 
from the LCDM simulation, while the densities are recovered with 
p-w relations from OCDM and TCDM simulations. 

2 on large scales, and even more on small scales. The effect of 
this model dependence could be reduced with the constraint 
on small scales, since it is possible to recover the small-scale 
power well by removing the saturations. 

Fig. [7| tests a different statistics - the one-point distri- 
bution function for the LCDM model. It is evident that the 
Gaussian profile leads to a drop of the probability at high 
densities. In other words, it statistically reduces the heights 
of density peaks. 



4 MAPPING THE POWER SPECTRA 

It seems that the inversion is no longer needed at least for 
determining the mass power spectrum if one establishes a di- 



10 



1 — ' — I — ' — I — ' 

□ Original "j 
• Recovered 




10 3 t" LCDM 
z = 3 

2 [_ Tj= 0.10 
J i I i L 



_L 







6 8 10 12 14 16 
P 



Figure 7. The one-point distribution function of the recovered 
densities (filled circles) and that of the original densities (open 
squares) . 



rect mapping betw een the flux and the mass power spectra 
l)Croft et al.l2002h . It is shown that given a set of cosmologi- 
cal parameters, there is a statistical mapping, which reliably 
recove rs the mass power spectrum from the flux power spec- 
trum ifOnedin fc Hamilt,onll20o'3l . even though the mapping 
is model dependent. 

The physical links between the three-dimensional linear 
mass power spectrum Vh(k) and the flux power spectrum 
can be summarised by the flowchart: 



P L (fc) — ■+ P NL (fc) 



^NL(k) 



P{k) ft(fc), 



where VNh(k) is the 3D non- linear mass power spectrum 
jPeacock fc Doddslll995 iTavlor fc Hamiltonlll996h. V§i.(k) 
is VNh(k) in redsh ift-space feaiserl Il987l : IPeacockl Il999l : 
IZhan fc Fanell2003lL and d 2 (k) = P F (k)/P(k). Notice that 
PnlCO i s anisotropic. We show in Fig.|H|that the non- linear 
evolution and the redshift distortion are significant at z — 3. 
The fact that the departure from linear evolution is very sim- 
ilar in both large and small box simulations indicates that 
such departure is real, and the cosmic density field has al- 
ready gone non-linear below 10 ft M pc at z = 3 (see also 
IZhan et al.ll 200 it IPichon et all 1200 J) . The angularly aver- 
aged 3D power spectrum V§s^{k) does not give a complete 
view of the difference between the real-space power spec- 
trum and redshift-space power spectrum, but it does show 
that peculiar velocities boost the power on large scales and 
reduce it substantially on small scal es. A two-dimensiona l 
projection of pN L (k) can be found in lPeacock et alJ i200ll) . 



In contrast, the mapping in ICroft et alJ J2002f) follows 
a simplified path: 



P L (fc) b JS P F (k), 



wher e Ph(k) = V L (y)ydy flKaiser fc Peacock! 

ll99llL and b 2 (k) = P F (k)/P h (k). Since V§ L (k) is 
anisotropic, P(k) and ^^(k) do not follow the same inte- 
gral. Given a one-dimensional density field in redshift space 
at z = 3, one can only measure P{k), and so Ph(k) is not ob- 
servable unless the density field is linear on all scales and no 
peculiar velocity is present. In other words, P(k) is directly 
connected to Pp(fc), while Ph{k) is not. 

Around the cosmic mean matter density, i.e. |<5| <g; 1, 
equation Q is approximately 
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Figure 8. The 3D power spectrum of the LCDM model at z = 3. 
The dashed lines are the real-space power spectra, the dotted 
lines are the redshift-space power spectra averaged over solid an- 
gles, and the thin solid line is the initial power spectrum (at 
z = 15) evolved to z = 3 by linear theory. The thin dashed line 
and the thin dotted line are from an additional simulation with 
the same parameters as the LCDM model except that the box is 
51.2 ft. - x Mpc each side. 

F(s) ~ e~ A - Af3e~ A 5(s), (7) 

where we have included the dependence on the redshift co- 
ordinate s explicitly. It is seen in the linearized equation 
10 that the Fourier modes of the flux are proportional to 
those of the overdensity when k 7^ 0. Hence, the flux power 
spectrum is proportional to the mass power spectrum 

P F (k) ~ A 2 f3 2 e- 2A P(k), k ^ 0. (8) 

Caution should be taken when using equation (J7J , because it 
is valid only if one smoothes the Lyq forest over large sc ales 
(see, for example. Izhan et aUfiofJll : IZhan fc Fandl2002ft . 

It is obvious from equation Q and Fig. Qthat the con- 
dition \S(s)\ <C 1 is never met where F ~ 1 or F ~ 0. There- 
fore, one has to include higher order terms of S(s). Taking 
the simplest case, in which S(s) = sin(fcos) and only 5 2 (s) 
term is added to equation (j), one immediately finds that 
Pp(fc) contains spurious power on the mode 2ko that P(k) 
does not contain. In general, the non-linear density-flux re- 
lation distorts, spreads, and mixes power in the cosmic den- 
sity field over different scales in the flux power spectrum. 
Other sources of dis tortion originate from linear filtering 
jGnedin fe Huilll998l) . line profile (see lGravllT99l3) . and in- 
strumentation. We do not focus on these effects, since they 
are well studied in the literature. 

Fig- |H]sliows the ratio b(k) and d(k). The dispersion in 
d(k) is larger than that of b(k), because both Pp(fc) and P(k) 
contribute to the dispersion of d(k) , while b(k) receives only 
a single contribution from Pp(k). Even so, the error prop- 
agation would cause a large scatter in determining P^Qt). 
In fact, the scatter can be large enough that one would not 
be able to determine the cosmological model based only on 
the recovered Ph(k). On the other hand, the model depen- 
dence of b(k) is strong enough, so that one may have to 
assume a priori the cosmological model to recover the mass 
power spectrum from the flux power spectrum. The differ- 
ence among the models is not solely due to a$. For exam- 
ple, the standard cold dark matter (SCDM) model has a 
slightly higher erg than the TCDM model, so one would ex- 
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Figure 9. The ratio b(k) = y/ ]'}.(k)/ P].(k) (upper panel) and 

d(k) = ^/P F (k)/P(k) (lower panel). The Lyo forests are divided 
into groups, each of which consists of 20 LOS samples. The power 
spectra are averaged within each group except Pi J (k), which is 
theoretical. The error bars of the LCDM model are la dispersions 
among 610 groups. Error bars of other models are comparable to 
that of the LCDM model. 

pect b(k) of the SCDM model to have the same shape as that 
of the TCDM model but with a shift. This is not observed in 
Fig. [5] however. It is also evident from the behaviour of b(k) 
and d(k) that the power spectrum of the underlying density 
field P(k) is significantly lower than the linear mass power 
spectrum Ph(k) on scales k > 0.02 (km s — ) . This is due 
mostly to the peculiar velocity, or redshift distortion. 



5 DISCUSSION AND CONCLUSIONS 

We have provided a scaling relation between the mean mat- 
ter density and the width of the saturated Lya absorptions. 
It has not been tested against variations of A or j3 [see equa- 
tion (2}]j because the p-w relation is the property of the den- 
sity field, while changing A or j3 can only affect the width 
of the saturation, and consequently determines which part 
of the scaling to be measured. This is also why the p-w re- 
lation slightly depends on the threshold flux, below which 
the underlying density cannot be directly extracted from the 
flux. 

The inversion with the p-w relation is able to re- 
cover the mass power spectrum fairly well on scales above 
2 /i _1 Mpc, but underestimates the power on small scales due 
to the use of the Gaussian profile. An improvement is sug- 
gested based on the observation that the power spectrum, 
with saturations removed, recovers the small-scale power 
well. Thus, by combining both the p-w inversion and the 
result after removing the saturations, one could get a good 
estimate of the power spectrum on all scales. 

The Gaussian profile does not preserve the high-density 
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tail of the one-point distribution function either. A more re- 
alistic density profile may improve the one-point distribution 
function, but the complication of substructures and peculiar 
velocities could result in only a small margin for improve- 
ment. 

The p-w relation provides an important constraint to 
the inversion of the Lya forest. To incorporate it in any 
inversion scheme, one needs to determine the p-w relation 
according to spectral resolution and noise. The threshold 
flux can be conveniently set to a small value above the noise 
level. That is why we have not added any noise in the sim- 
ulations: noise outside the saturation can be routinely dealt 
with, while noise in the saturation, being below the thresh- 
old, does not matter. 

The density is so far recovered in redshift space. 
The real-space density can be obtained iteratively 
iNusser fc Haehnelli Il999t) . and the real-space three- 
dimensional ma ss power spectrum can be recovered using 
matr ix method ll Hui|| 19991) . or less rigorously by differentia- 
tion JCroft et al.lll998ft . 

It is discussed in Section ^. ll that the pseudo-hydro tech- 
nique is reliable above 200 h~ kpc at z = 3. This means 
that while the bulk of our analysis is valid, the details below 
200 h~ kpc may be subject to hydrodynamic effects. In the 
future, we will put the p-w relation on test using hydrody- 
namical simulations. 

Finally we find that the non-linear evolution and red- 
shift distortion have significantly changed the mass power 
spectrum at z = 3. Therefore, the linear mass power spec- 
trum Pl (k) is not directly observable, and it is substantially 
larger than the mass power spectrum P(k) on scales k > 0.02 
(km s^ 1 ) -1 at z = 3. 
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